4Ф 


УДК 681.3 


Л.П. Фельдман 
Донецкий национальный технический университет, Украина 
Ге|атап(@г5.4оа.ЧопеК ма 


Устойчивые одношаговые блочные методы 
численного решения жестких обыкновенных 
дифференциальных уравнений 


В статье предложены устойчивые одношаговые блочные методы решения задачи Коши для систем 
обыкновенных дифференциальных уравнений. Методы обладают высокими характеристиками устойчивости, 
что позволяет применять их для решения жестких задач, и легко отображаются на параллельные 
структуры произвольной архитектуры. 


Статья содержит обобщение результатов исследований [1-5], посвященных ус- 
тойчивости традиционных последовательных методов численного решения задачи 
Коши жестких систем обыкновенных дифференциальных уравнений и является про- 
должением ранее опубликованных работ [6-9]. В ней рассматривается устойчивость 
параллельных численных методов решения задачи Коши жестких систем обыкновенных 
дифференциальных уравнений. Предлагаемые блочные одношаговые многоточечные 
разностные схемы обладают высокой точностью, легко распараллеливаются и жестко 
устойчивы. 


Постановка задачи 


1. Рассматривается устойчивость решения задачи Коши жестких систем обык- 
новенных дифференциальных уравнений 


4х 


— = Л, х) >, х( 1) =х, (1) 
Ге 
одношаговым многоточечным методом, определяемым формулой: 
К — И 
Ил = Ин, 0+ 0+ У а, Ль; 51ЕЪК, ПЕЬМ. (2) 
7=1 


При исследовании устойчивости блочных разностных методов для жестких систем 
уравнений, так же, как и для классических методов, обычно рассматривают модель- 
ное уравнение (3) 

ах — 

@ 
где /, — произвольное комплексное число. Свойства различных методов анализиру- 
ют на примере модельного уравнения (3). 


Ах, (3) 
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Для модельного уравнения (3) разностные уравнения (2) примут вид 


[4 


гаи о +1 Ви, 0+ Уаз уи,; |) 1=ЪК, п=ЬМ. (4) 


НЕ 


Обозначим матрицу коэффициентов и вектор-столбец системы уравнений (2) через 


Аа. ВФ ЕЕ, 


Представим систему уравнений (3) в векторной форме. Для этого введем 


векторы 
= ше аа 


п 


В векторной форме уравнение (3) будет иметь вид 
(Е- ии, =и„о(е + иБ), и = Ат. 
Определив из (5) И» ‚ получим 
И, =(Е- и) (е + иБ)иь=5и, в, 
где вектор 5 определяется выражением 


5=(Е- мА) (е +05), 


ИЛИ 
5 
1- дал, —иа, >, На... Иа 1+ 46, 
х —иа,, 1- а, ›,- Ма» з,...— Мау 1+ дБ, 


(5) 
(6) 
(7) 
51 
52 
(8) 
5% 


Определим матрицу перехода от блока значений 0,_ = (и, 11.И,1,›..„, ‚К блоку 


О, „1. Для этого преобразуем полученную систему вида (6) к эквивалентной системе 


0,0,.....0,5) 
0,0,.....0,5) 

5= КИ, =50 ьт=Ь,... 
0,0,.....0,5 


(9) 


Таким образом, устойчивость численного метода (2) определяется собственными 
значениями матрицы 5 [2]. Видно, что элементы этой матрицы являются рациональ- 
ными функциями от д и собственные значения 4 матрицы 5 также будут зависеть 


от д. Поэтому необходимо найти те области для собственных значений д, в кото- 


рых (и) <. 


По виду матрицы 5 можно заключить, что она имеет А — / нулевых собственных 


чисел 4,/=4, =... =4,, =0и 4, =5,. Поскольку ранг матрицы равен единице, 


каждому собственному значению соответствуют различные ненулевые собственные 
векторы. Таким образом, функцией устойчивости метода (1) является 


Чк(И)=5ь- 


(10) 
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Из (7) следует, что 


1- дат, а дара 
Рей А] д - 2иа,,1- а>»,- Иа... 1+ ИБ 


РЕ а очень в 


5 


Матрица 5 имеет четыре ненулевых собственных вектора. Найдем область дейст- 
вительных значений д, для которых |4 [и] | =. 


60 + 120 т+ 105 п® + 50 мВ + 12 пи 
60 - 120 т+ 105 м - 50 п + 12 пи 


Ч4 = 


В последнем выражении модули числителя и знаменателя равны, следовательно, 
на мнимой оси Ч 4[1&] | = /. Поэтому областью устойчивости метода является вся левая 


полуплоскость, включая и мнимую ось. Таким образом, одношаговый четырехточеч- 
ный блочный метод А-устойчив. 
Проще определять функцию устойчивости, используя формулу (11). 
Определитель матрицы Е — ИА будет равен: 


РеНЕ-иА]=1-4и+7.58333/? —9 р? +7.4236 111 —4.45° +1.95311б —0.603968 и? +01. 


Заменим в матрице Е- А последний столбец вектором 5, получим матрицу А». Най- 


дем ее определитель 
Бе[А,]=1+4и+7.58333 и” +9 +7.4236 1 +4.45 и? +1.953114° + 0.603968 и” + 0.11111. 
Таким образом, функция устойчивости метода будет равна 


_ 1+41+7.58333м° +9” +7.42361* +4.45 и° +1.953 11° +0.603968 и” + 0.11111" 
1-4и+7.58333 и? -9 и? +7.42361и* —4.45 и +1.95311° —0.603968 и’ +0.11111* ` 


Ч [и] 


В рассмотренном примере модули числителя и знаменателя равны, следовательно, 
на мнимой оси | 8[1а] | = /. Поэтому областью устойчивости метода является вся левая 


полуплоскость, включая и мнимую ось. Таким образом, одношаговый восьмиточеч- 
ный блочный метод также А-устойчив. Очевидно, что одношаговые многоточечные 
блочные методы более высоких порядков также будут А-устойчивыми. 

Рассмотрим пример использования одношаговых многоточечных блочных ме- 
тодов. Типичным примером жесткой задачи является уравнение Ван-дер-Поля. Его 
можно записать в виде системы: 


ея 


жало (х? +) (12) 


Для расчетов положим @ = 100, х =2,х› =0. 
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Построим графики полученных решений: 
х2 


Х1 


Рисунок 1 — Фазовая траектория (предельный цикл) решения задачи 


Полученное решение свидетельствует о возможности решения жестких задач 
блочными одношаговыми многоточечными методами. На рис. 1 видна переходная 
зона от начального положения к предельному циклу. Ее так же можно получить по 
приведенной программе. На рис. 2 приведена переходная фаза, полученная с мень- 
шим шагом интегрирования. 


х1 


Рисунок 2 — Переходная фаза решения задачи (12) 


Заключение 


На основе проведенного анализа параллельных разностных методов можно сде- 
лать вывод, что одношаговые многоточечные разностные (1) методы являются А- 
устойчивыми. Таким образом, эти методы могут использоваться при численном 
решении жестких систем обыкновенных дифференциальных уравнений. Предложенная 
методика может быть использована при исследовании устойчивости параллельных 
одношаговых методов с любым числом точек в блоке. Проведенные численные решения 
одношаговыми блочными методами тестовых жестких систем практически подтвердили 
их надежность и эффективность. 
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Стйк! однокроков! блоков! методи чисельного розв?язання жорстких звичайних 
диференщальних рвнянь 

У стати запропоновано стйк! однокроков! блоков! методи розв’язання задач! Кош! для систем звичайних 
диференщальних равнянь. Методи мають висок! характеристики стйкост!, що дозволяе розв’язувати 
жорстк! задач! та легко в:дображаються на паралельн! структури дов1льно! архтектури. 


Г.Р. Е@4тап 
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ЗаЫе опе-5ер Моск те#о4$ оЁ питега| 4ес151оп оЁ СаисВу’$ ргоМет г зузет$ оЁ ог4тагу Аегепна1 
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